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I. INTRODUCTION 

The second law stipulates that the entropy of an iso- 
lated system cannot decrease. In a recent development, 
[l|, see also 0-0: it was suggested that the second law 
can in fact be split in two. The total entropy production 
(EP) Stot is the sum of two constitutive parts, namely 
a so-called adiabatic Sa and nonadiabatic Sna contribu- 
tion. These separate contributions arise from the fact 
that there are two mechanisms that lead to the time- 
symmetry breaking characteristic of a dissipative pro- 
cess, namely the application of steady nonequilibrium 
constraints (adiabatic contribution) or the presence of 
driving (nonadiabatic contribution). The crucial point 
is to note that each of these contributions is separately 
non-negative. We can thus identify "three faces" to the 
second law, the positive rate of production of the total 
EP, of the adiabatic EP and of the nonadiabatic EP: 



Stot>0 , ^„a>0 , Sa>0. 



(1) 



Explicit expressions for Stot, Sa and Sna, including a de- 
tailed mathematical and physical discussion, were given 
at the level of a master equation description in the pre- 
ceding companion paper [6[. However, in many applica- 
tions, a description on the basis of a Langevin or Fokker 
Planck equation is more appropriate. The purpose of this 
paper is to provide a detailed discussion of the adiabatic 
Sa and nonadiabatic Sna EP for such a description. It 
should be noted that the master equation description is 
the more general one, including the Langevin and Fokker- 
Planck case as a special limit. The rather technical tran- 
sition between both descriptions based on such a limiting 
procedure is given in the appendix of this paper. In order 
to be both self-contained and physically motivated, we 
derive the same expressions for the adiabatic and nonadi- 
abatic EP directly from the Fokker-Planck equation itself 
in the main text. The application to Langevin equations 
driven by Gaussian white noise is immediate since there is 
a mathematical equivalence between Fokker-Planck and 
Langevin descriptions [7] . Note finally that we focus here 
on the EP rates for which we provide explicit expres- 
sions. We have previously obtained the results for the 



time-integrated trajectory-dependent versions of the var- 
ious EP contributions [l|. The latter can be expressed 
in terms of relative entropies between probabilities for 
paths in a direct and various types of reverse experiments. 
While these results have a profound meaning by reveal- 
ing the temporal symmetry breaking associated to each 
contribution, the expressions for the EP rates given here 
do not refer to any reverse experiment and are thus much 
easier to measure or calculate. 



II. FOKKER PLANCK EQUATION 

A. Total entropy balance 

Our starting point is the Fokker-Planck equation de- 
scribing the time evolution of the probability density 
Pt — Ptix) for the variables x: 



Pt = -dxJt 



with 



j^-T^Ji 



H 



utpt - DtdxPt 



M 



M 



jr^urpt-nrd^pt 
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Ut 
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(3) 

(4) 
(5) 



The quantities pt, Jt, Jt ,''^t , ^t , 6tc. are all func- 
tions of the state x of the system, although this is not 
written explicitly for simplicity of notation. The proba- 
bility density changes in time due to different processes 
v. As a result a probability fiux J^' is associated to 
each process. We also assume general time-dependent 
drift and diffusion coefficients Mj and D^ . 

We now proceed with the identification of various EP 
contributions, whose positivity can be guaranteed on a 
purely mathematical basis. As we will see in the next 
section, these terms can be interpreted as genuine forms 
of EP when the proper physical content is taken into 



account. We focus on the time evolution of the Shannon 
entropy for the system: 



S{t) = - dx pt \npt 



(6) 



Note that we assume here and henceforth that the Boltz- 
mann constant ks is set equal to 1. Using ([2]) one finds: 



S{t) = - dx Pt Inpt 



Y.hxj[^^'-^ 



J2jdxJ, 



ri^) 



Ut 






(7) 



We used partial integration (assuming that the boundary 
contributions cancel) and the definition of the probability 
flux ((4]) to get the results on the second and third line. 
One identifies the following two parts: 



S[t)^Se{t)+S^{t). 



(8) 



The first term corresponds to the entropy flow into the 
system: 



Se{t) 



^I'-'''W 



(9) 



The second term is the non-negative irreversible EP: 



{Jtr 






(10) 



These expressions are in agreement with the results given 
previously in the literature for the case of a single process 

We mention two additional properties of the total EP. 
First, we note the following inequality: 



iJV 
S^{t) > I dx^^ > 0, 



(11) 



which shows that the total EP is underestimated if the 
constituent processes v are not properly identified. The 
above result follows from the following inequality, valid 
for any set of numbers j/i > 0: 



Y^ ^ ^ {z2i^i 






Vi 



J2^y^ 



(12) 



To proof this inequality, consider first the case with all 
Xi > 0. The above inequality is identical to Jensen's in- 
equality {x/y){y/x) > 1, where the averages are over the 
variables Xi/yi and yt/xi with respect to the probability 
distribution pi ~ Xif^^Xi. The above inequality will 
hold a fortiori if not all Xi > 0, since the l.h.s. is insen- 
sitive to a change of sign of the variables a;,, while the 
r.h.s. can only become smaller. 



Second, the irreversible total EP has the familiar form 
of a sum over fluxes times forces: 



^,(t) = ^|dxjf)x('^), 



with the force associated to process v given by: 



Dr{x)pt{x) 



(13) 



(14) 



B. Thermodynamic interpretation 

The above description can be postulated on a purely 
phenomenological or mathematical basis, as it corre- 
sponds to the general equation of evolution for a con- 
tinuous Markov process. We now make a number of 
comments that validate the model and the derived ex- 
pressions for the various types of EP, from a physical 
point of view. 

A Markovian stochastic evolution for the degree of free- 
dom of interest, the system, originates from the elimina- 
tion of fast degrees of freedom, the reservoir, that do not 
need to be described because they are at instantaneous 
equilibrium for any given state of the relevant variables. 
Furthermore, it is assumed that the different processes 
u corresponding to groups of fast eliminated variables 
(reservoirs) at different equilibrium values, do not di- 
rectly interfere which each other. Otherwise this interac- 
tion would need to be described and the set of variables 
describing completely the physical process would need to 
be enlarged. The correct identification of the reservoirs 
I' is a crucial step since we have seen in (jlip that fail- 
ure to do so will typically lead to an underestimation of 
the EP. At constant (in time) drift and diffusion coeffi- 
cient, when all reservoirs but one are switched off, say 
V, the system should reach an equilibrium steady state 
distribution corresponding to the thermodynamic prop- 
erties of the reservoir v and satisfying the condition of 
detailed balance X^ = 0. In presence of different reser- 
voirs V, the steady state is out of equilibrium because 
it will break detailed balance. Indeed, all the reservoirs 
try unsuccessfully to impose their equilibrium value on 
the system. The essential step to connect the stochastic 
description to the present thermodynamics discussion is 
the local detail balance condition: 



M) 



D 



M 



P^^^Ft. 



(15) 



This relation implies to identify the energy Et of the sys- 
tem, because the force Ft is the negative derivative of 
the energy: Ft = —dxEt (for more details, see also ap- 
pendix and Eq.(22) in [6] for the corresponding relation 
for the master equation). It also mathematically guar- 
anties that in presence of a single reservoir v and of a 
time-independent external force F , the probability dis- 
tribution of the Fokker-Plank equation ^ will eventu- 
ally reach the equilibrium distribution p'^'^ ^ exp {—^E}. 



Since p5)) translates the fact that each of the reservoirs 
remains at equiUbrium, the EP P^ is also the total EP 
since no irreversible processes take place in the reservoirs: 
Si{t) = Stot{t)- For the same reason, the entropy flow (jH]) 
into the system corresponds to minus the entropy change 
into the reservoirs: Se{t) = —Sr{t). 

By introducing the generalized mobility /ij of process 



i^ as Mj = Ht Ft, we see that p5|) is in fact the gen- 
eralized fluctuation-dissipation Einstein relation fj,^ = 
fjMjJY' . This clarifies the thermodynamic meaning of 
the force ([Tl)) since J^' /pt can be identified as the local 
speed wt"^' of process v so that x[''^ = z;('')/(/^('')t('')) 
corresponds to the usual thermodynamic expression of a 
force, speed divided by mobility, over the temperature. 
Furthermore, using (fT5|) in ([9]), we find that the entropy 
flow takes the familiar thermodynamic form 

where heat flowing into the system is given by 



Q^'Kt) 



dx A^^Ft 



Introducing the system energy 

E{t)= [ dxptEt, 



(17) 



(18) 



we find (using integration by part and neglecting the 
boundary terms) that the first law of thermodynamics 
assumes the familiar form: 



E{t) = w{t) + Y,Q(''\t), 

where the work is given by 



W{t) 



dx ptEt 



(19) 



(20) 



We have thus shown that the local detail balance 
condition (jlSp provides an explicit connection to ther- 
modynamics and justifies the names used for the various 
entropies in the previous section. 

We should note however that for systems subjected to 
nonconservative forces, the local detail balance condition 
(ITSI) will not be satisfied. Even in presence of a single 
reservoir v, the steady state will break detailed balance 
and will thus be a nonequilibrium steady state. 



C. Adiabatic and nonadiabatic entropy balance 

To identify a component related to the relaxation of 
the system, we introduce the instantaneous steady state 
solution pf, being the normalized (supposedly unique) 
solution of the following equation: 



da^Jf 



0, 



(21) 



with 



st(ii) 



J'' = E J' 



"'' DApf 



't - / ^'H - UtPt 

jr'=uY'pf-D't'^d,pf. (22) 



It corresponds to the the steady state solution of the 
Fokker Planck equation ([2]) , if the drift and diffusion co- 
efficients are frozen at their instantaneous values. We 
can now rewrite the fiux Jj , cf. (|3|), as follows: 

j{v) jst{u) 



Pt 



pf 



-D['^dA\n^ 



Pt 

pf 



(23) 



As a result, the following expression is identically zero: 



E / ^^-^t 



st(u) Ph 
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dxjfd^ ( 4 
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Pt 



(24) 



The last step in ((24|) follows from (|21|) by partial inte- 
gration, assuming that the boundary term is zero (infi- 
nite system or system with periodic boundary condition) . 
The separation of total EP (flUl) in two contributions that 
are separately positive is now straightforward. We write 
the integrand in ([TU]) as follows: 
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(25) 



From (1211) and (gHI), we conclude that the total EP CT 
can be written as follows: 



Stotit)=Snait)+Sait), 



(26) 



with the following explicit expressions for the nonadia- 
batic and adiabatic rate of EP: 



Sna{t) ^y^ldx 



/ t(!^) jst(v)\ 2 

Pt I 'Jt '^t 



Di"^ \ Pt 



pf 



Pt f Ji 



D 



jst(v) \ 2 



pf 



>0 



> (27) 
(28) 



These quantities are clearly non-negative. The nonadi- 
abatic EP is zero for an infinitely fast relaxing system 
being all the time in the instantaneous steady state. The 
expression for the adiabatic EP is similar to the total EP, 
but with the steady state contributions J^ /pf rather 

than the actual jf /pt singled out. Both expressions 
can also be obtained as the limits of the corresponding 
expressions for the nonadiabatic and adiabatic EP for a 
master equation, cf. appendix [Xj In case of a system in 
contact with a single reservoir and subjected to a noncon- 
servative force (the steady state break detailed balance), 
the adiabatic EP is the housekeeping heat divided by the 
reservoir temperature 0, 0: I12h14| . 



We now make a number of further comments includ- 
ing alternative expressions for the adiabatic and nonadi- 
abatic EP. First, from 



Pt 



1 

Dt 


yPt 


JSt 

•^t 
Pf 


J 


1 


(Jt" 

' V Pt 


) 


jst(v) 

•Jt 


Di"' 


pf 



= -dJin^] (29) 



it follows that the nonadiabatic EP can be written in 
terms of compound quantities (obtained by summation 
over the processes v) only: 



^na\t) 



dx 



Pt fJt J! 



t\ 2 



A \pt pf 
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(30) 



We conclude that the nonadiabatic EP is not sensitive 
to the identification of the various separate processes 
describing the exchange with different reservoirs. The 
rationale is that this EP reflects relaxation EP within 
the system itself. We recall furthermore that the total 
EP, which is the sum of the adiabatic and nonadiabatic 
contribution, can only decrease upon coarse-graining the 
processes i'. Hence, we conclude that the adiabatic EP 
is underestimated in exactly the same way as the total 
EP, when the constituting processes v are not properly 
identified. 

Second, we mention the following alternative form of 
the nonadiabatic EP: 



^na(i) — 



dxptln—T; 
pf 



(31) 



obtained from (pO)) with ((29)) . This form is quite con- 
venient for a direct calculation of the nonadiabatic EP 
when the probability distribution pt is known explicitly. 
We also note that for constant in time drift and diffu- 
sion coefficients, one has that H = —Sna < 0, where 
H = J dx Pt Inpt/pf > is a Lyapunov function. This 
thus proves the convergence of pt to the (supposedly 
unique) steady state p** |15[ . 

Third, we introduce another entropic contribution, the 
so-called excess heat UM: 



^ex\t) 



dx Jfd, 



Inpf. 



(32) 



This expression allows to complement the familiar EP 
balance equation ([8|) with two other balance equations, 
leading to an alternative presentation of the "three faces" 
of the second law. Indeed, one immediately verifies that 



S{t) = -S,At)+S^a{t) 



(33) 

(34) 



Each of these balance equations features the sum of an 
exchange term, the excess entropy, with has no definite 
sign, plus an irreversible non-negative EP term. The 
nonadiabatic term is related to the system properties and 
is independent from the constituting processes ly. This 



is not the case of the adiabatic term which represents 
the dissipation incurred via the contacts with the various 
reservoirs. We note that when considering transitions 
between steady states, ([55| becomes the second law of 
steady state thermodynamics [l^, \lM ■ 

Fourth, upon introducing the following thermody- 
namic forces: 



X^t"^ = ^ - 4"^ + Nt 



Dl-"pt 



(35) 
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each of the irreversible EP terms can be written under 
the familiar form of a sum over fluxes times forces: 



(i') /.(") 



Snait) = / dx JtNi 



(37) 
(38) 



III. APPLICATIONS 

A. Brownian particle in contact with two thermal 
reservoirs 

We consider an underdamped Brownian particle in 
contact with two separate heat baths at temperature T^^^ 
and T^^) . Such a model has been studied in the context of 
an analysis of the Feynmann ratchet [la ] . It corresponds 
to the simplest model for thermal conduction by a single 
degree of freedom. The more suggestive presentation is 
via an equation of motion written under the form of a 
Langevin equation: 



-(7^^ 



■It 



(2) 



)v + ^/2 



2jI'^tW^W + A/27f ¥(2)^(2) 



(39) 



with ^(^^ and ^'-^^ independent Gaussian white noises of 
intensity 1 . The variable v plays the role of the " speed" 
of the Brownian particle (although such a speed does not 
exist in a mathematical sense, only its increments are 
well defined), while 7} and 7t are the friction coeffi- 
cients appearing due to the contact with the respective 
reservoirs 1 and 2 which we assume externally control- 
lable (even if this might be physically not very realistic 
it serves to illustatc our results) . The mass of the particle 
is taken equal to unity. Note that we do not take into 
account a spatial degree of freedom. This further simpli- 
fication corresponds to a thermal contact tightly bound 
to a specific location (in a "delta function" potential). 

In the context of the Fokker Planck description, we 
identify the following two drift and diffusion coefficients 
(i. = l,2): 






-It v 



D 



H 



It 
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(40) 



Note that we have incorporated the appropriate fluctua- 
tion dissipation theorem through the relation linking the 
diffusion to the friction coefficient (here, E = v"^ /2 cor- 
responds to the kinetic energy). We mention a further 
peculiarity of this model. The Langevin equation can be 
rewritten as: 



with 



ut = 



E 



-Jtv 
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a^Ea''^^-^*^*- 



(41) 



(42) 



Consequently, if the distinction between the two pro- 
cesses v is not made, this situation corresponds to Brow- 
nian particle in contact with a single heat bath at tem- 
perature 

^* - (1)^ (2) (43) 

and the steady state corresponds to equilibrium (zero EP, 
with equipartition {v'^)t — Tt, see also below). This point 
illustrates our discussion concerning the physical input 
needed to validate the expression for EP, and in partic- 
ular the correct identification of the basic processes that 
are taking place. As mentioned before, the nonadiabatic 
EP will be correctly reproduced, but both total and adi- 
abatic EPs will be underestimated. In the present case, 
the underestimation is dramatic since the stationary dis- 
tribution of the reduced description corresponds to ther- 
mal equilibrium, so that the coarse-grained adiabatic EP 
will be identically zero. 

The probability distribution pi(w) for the speed v obeys 
the following Fokker-Plank equation: 



Pt(w) = ltdv[vpt{v) + Ttd-uPtiv)]. 



(44) 



The solution to this equation is a Gaussian distribution, 
if it is so initially: 



Pt{v) 
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yfhd^^ 



■ exp 



2(w2 



(45) 



Its time evolution is completely determined by that of 
the first and second moments, obeying the following set 
of equations: 

dt{v)t = -lt{v)t 
dt{v^)t^~2jti{v^)t-Tt). 

The steady state form is given by 

1 / V 
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One easily verifies that: 
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It is also convenient to define Q{t) = Q^'^\t) + Q^'^\t). 
Inserting the above results in the expression for the var- 
ious forms of irreversible EP, we find: 



Sa{t) 



Q{t) 

Tt 



E- 






Sna{t)^Q{t){- 



v')t 



i) 



7 Tt 
lT{v^)t 



T1T2 



Stotit) ^ Sait) + Snait) 



(50) 



For the initial confition pt=oiv) = S{v) and in absence of 
externally control of the friction coefficients, the equation 
for the second moment reads: 



{v')t 



T(l 



-27^ 



(51) 



which lead to the simplification: 



Sait) = ^^^il-e- 
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' TIT2 



-47* 



^«a(0= 7^—^-27* 
Stotit) ^ Sa{t) + Sna{t) 



(52) 



These EPs reproduce expected properties. In absence of 
external driving the nonadiabatic EP decays to zero as 
the steady state (where (w^) = T) is approached. The 
adiabatic contribution associated to the application of 
the non-equilibrium boundary conditions tends toward 
the usual thermodynamic expression for the EP associ- 
ated to a steady heat flux between two reservoirs by a 
device with a thermal conductivity equal to k = 7172/7 
|16| . We finally note that Fourier law is recovered in the 
steady state Qi = —Q2 = n{Ti — T2). 



B. Driven Brownian particle on a circle 



We next consider an overdamped Brownian particle 



x = ut + V2D^, 



(53) 



where ut and D are the (time dependent) drift and 
(time independent) diffusion coefficient both position in- 
dependent. We furthermore assume x € [0, 1] with peri- 
odic boundary conditions. This model was discussed in 
stochastic thermodynamics as a simple example of a par- 
ticle driven by a nonconservative force 17]. For simplic- 
ity, we consider the special initial condition pt=Q — S{x). 
The exact time-dependent solution of the Fokker Planck 
equation is expressed in terms of the well-known solution 
on the infinite line: 



namely (x G [0, 1]): 
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(55) 



It converges to the steady state solution: 
pf{x) = l. 



(56) 



This model has a peculiarity: the steady state distribu- 
tion is identical to the equilibrium distribution pf'{x) = 
p^'^{x) — 1. As a consequence one has that J|* = Ut. Fur- 
thermore, once pt has relaxed to this distribution, this 
remains the case even while a time-dependent driving ut 
is still applied. We find that the various EPs read: 



Sa{t) ^ ^ 



Stotit) ^ Sa{t) + S„a{t) 
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The adiabatic EP is proportional to the square of the 
externally applied drift. The nonadiabatic EP is given by 
a more complicated expression, but goes to zero as the 
probability distribution relaxes to the uniform stationary 
distribution. It remains zero once this distribution is 
reached, even if a time-dependent driving ut persists. As 
an illustration, we reproduce the results for the various 
EP contributions in Fig. ([1]), including the effect of a 
switch in the driving speed. 



IV. CONCLUSION 




In this paper, we have identified the non-negative EP 
as well as its two non-zero contributions, the adiabatic 
and the nonadiabatic part, for Fokker-Planck dynamics. 
This parallels a similar identification for master equa- 
tion dynamics presented in the companion paper [6| . We 
have shown that this identification allows to "split the 
second law in two parts" . It remains to be seen what are 
the implications of this "doubling" of the second law. In 
particular, we speculate that it should imply the impos- 
sibility of some physical phenomena, being incompatible 
with the inequalities, that it may provide novel limits, for 
example on efficiencies of machines, or may be linked to 
novel symmetries, such as the symmetry of Onsager co- 
efficients. Finally, we reiterate that the positivity of the 
adiabatic and the nonadiabatic EP obtained here follows 
by Jensen's inequality from the detailed fluctuation the- 
orems derived in [1[ . The latter deal with the trajectory- 
dependent adiabatic and nonadiabatic EP, and thus re- 
veal a much more detailed and deeper statistical symme- 
try deriving from micro-reversibility. 



Appendix A: Fokker-Plank limit of the master 
equation 

In the companion paper [6|, we derive the adiabatic 
and nonadiabatic EP for Markov processes obeying a 
master equation (cf. eq. (7)). The results for the Fokker 
Planck equation given in the main text can be derived 
by applying an appropriate limiting procedure, similar 
to that of Ref. [18|. As a starting point it suffices to 



FIG. 1. (Color online) Probability distribution at difTerent 
times and total, nonadiabatic and adiabatic EP {D = 0.02). 
The initial divergence of the nonadiabatic EP is due to the 
singular initial condition pt^o = 5{x). As the uniform dis- 
tribution is approached, the nonadiabatic EP decreases. The 
application of a switch from the initial value 0.8 to the value 
1.2 of the drift, cf. inset, has no effect on the nonadiabatic 
EP, but results in an additional adiabatic EP. 



consider the case of a tri-diagonal transition matrix, i.e., 
the only non-zero non-diagonal elements of W^^i mi^t) 
are those with m' — m ± 1. The master equation thus 
has the following form: 

p- = - E (•^ili.mW - 4r,L-iW), (Ai) 



where 
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The idea is that the (general) nearest neighbor random 
walk in the variable m goes over into a (general) dif- 
fusion process for a continous variable x — me. We 
illustrate the procedure for x G [—L, L] with reflect- 
ing boundary conditions, covering in the limit L — >■ oo 
the case of real variables. A similar procedure can be 



applied for perfodic boundary conditfons. We consider 
m = -N,-{N - l),---0, 1,2,--- ,A^, with reflecting 
boundary conditions, Wn+i,n = W-n-i,-n = 0. We 
take the limit e — > and N — >■ oo, where N — L/e, with L 
fixed, obtaining a continuous variable x = me € [—L,L]. 
Using exp{±drn}fm = /mil and 



Pt =p{x,t) ^p^{t)/e 



(A4) 
(A5) 
(A6) 
(A7) 



we find that (|A1[) goes over into the Fokker-Plank equa- 
tion: 
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Y.d.J^ 
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(A8) 



Using these results, it is easy to verify that the various 
EP from section II of Q lead to the EP of section II of 
the present paper. 



Finally we note that the local detailed balance condi- 
tion with respect to the various processes v given by eq. 
(22) in [6| reduces in the Fokker-Plank limit to (fTS)) . 
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with zero flux boundary conditions. Similarly the forces 
(28)-(30) in fsj become 
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(All) 
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